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1. Introduction 


Deep in the deconfmed phase of QCD, with either temperature T or quark chemical potential 
jti much larger than the QCD scale Aqcd, asymptotic freedom eventually leads to a strong coupling 
constant g that is sufficiently small to permit the use of perturbation theory for calculating the 
thermodynamics of quarks and gluons under extreme conditions. Indeed, much effort has been 
invested in calculating the thermodynamic potential of QCD at high temperature. It is by now 
known up to and including order g 6 log(g) mm- These calculations require resummations 
of ordinary perturbation theory to cope with infrared problems, which is done most elegantly and 
efficiently by means of effective field theories [|5]]. Up to and including order g 6 log(g) one can 
in fact set up an expansion into Feynman diagrams, but this breaks down at the order g 6 , which 
receives nonperturbative contributions from the chromomagnetostatic sector. However, already 
the expansion up to that point suffers from unusually poor convergence properties and a strong 
dependence on the renormalization point even at temperatures many orders of magnitude higher 
than Aqcd> which seems to reduce calculations of this kind to a purely academic enterprise. The 
perturbative expansion appears to become reliable only for coupling constants so small that the 
thermodynamic potential is anyway very well approximated by the Stefan-Boltzmann result for a 
noninteracting plasma. 

However, this problem is not specific to QCD at high temperature with its nonperturbative 
magnetostatic sector. Similarly poor convergence behaviour appears also in simple scalar field 
theory [Q], and even in the case of large-A <j) 4 theory [7]|, where all interactions can be resummed in 
a local thermal mass term. As soon as one starts expanding out in a series of powers and logarithms 
of the coupling, the result for the thermodynamic potential goes wild, whereas the exact result that 
is available for the large-A <j) 4 theory is completely unspectacular and smooth. 

Various techniques have been developed to restore the convergence of the perturbative series. 
The more physically (rather than mathematically) motivated ones include “screened perturbation 
theory” [X. l f|. its generalization to gauge theories (“HTL perturbation theory”) [[To, 11, f2|. [i~3| |. and 
the use of 2PI techniques which put the emphasis on a quasiparticle description [|l4|, 15, l6|, 17]. In 
particular the latter results agree very well with available lattice data down to temperatures about 
three times the phase transition temperature so that the range of applicability of weak-coupling re¬ 
sults does not appear to be restricted to uninterestingly high temperatures. Moreover, these results 
also agree well in this regime with the perturbative results of Ref. []3|, j4|] to order g 5 or even g 6 log(g) 
provided the effective-field-theory parameters used in this approach are not treated strictly pertur- 
batively, but are kept in the form in which they appear naturally jl8|]. The difference this makes 
is shown in Fig. [I], where the light-gray band shows the renormalization-scale dependence (a mea¬ 
sure of the theoretical uncertainty) to strictly order g 5 , and the medium-gray band shows that of the 
unexpanded three-loop result, both compared to the lattice result of Ref. & represented by the 
thick dark-gray curve. In the unexpanded three-loop result, the renormalization scale dependence 
is highly nonlinear such that an extremum is attained at the upper boundary marked PMS for prin¬ 
cipal of minimal sensitivity 1 ; a different optimization of the perturbative result is fastest apparent 
convergence (FAC), which turns out to be close to the former. 


1 A slightly different implementation of the PMS to the thermodynamical potential with roughly comparable results 
has recently been presented in Ref. (§]■ 


2 


















Thermodynamics of QCD at large quark chemical potential 


Anton Rebhan 



T/T c 


Figure 1: Three-loop pressure in pure-glue QCD normalized to the ideal-gas value with unexpanded 
effective-field-theory parameters when ft is varied between nT and AnT (medium-gray band). The broad 
light-gray band underneath is the strictly perturbative result to order g 5 with the same scale variations. The 
full line gives the result upon extremalization (PMS) with respect to ll (which does not have solutions be¬ 
low ~ 1.37)); the dash-dotted line corresponds to fastest apparent convergence (FAC) in m 
jl ss \ .19nT. (From [ 181) 


which sets 


In the following, we shall first consider the extension to nonzero quark chemical potential 
of the perturbative results obtained by Vuorinen by means of the effective field theory provided 
by the technique of dimensional reduction. After a comparison with lattice results on the one 
hand, and analytical results from the large-tVy limit of QCD on the other hand, we shall turn to the 
case of chemical potentials much larger than the temperature. This regime turns out to be much 
richer than the one at high temperature. For sufficiently low temperature, new phases with colour 
superconductivity occur which can in fact be analysed by weak-coupling techniques applied to full 
QCD [21], 22, 23, 24j 25, 26, 27||. In this regime, and also in the normal phase at temperatures above 
the critical one, non-Fermi-liquid effects play a crucial role. As we shall describe, the relevant 
effective field theory in this situation is those of hard dense loops (HDL) [ ]28] . [29j], and we have 
used it to calculate systematically non-Fermi-liquid effects in entropy and specific heat as well 
as in the dispersion laws of fermionic quasiparticles at low temperature and high quark chemical 
potential. 


2. Dimensional reduction at high temperature 

When the coupling constant is small, g <C 1, and temperature T is the largest scale, there is 

a natural hierarchy of scales: the “hard” scale T, the “soft” scale gT, “ultrasoft” g 2 T, _In the 

imaginary-time formalism, the hard scale T is carried by nonzero Matsubara frequencies co n = 
ninT with n even/odd for bosons/fermions, and all the softer scales involve only the zero modes 
of the bosons. The “soft” scale gT is where collective phenomena such as Debye screening and 
Landau damping occur 2 and its effective field theory is Yang-Mills theory reduced to three spatial 

2 In plasmas with a momentum-space anisotropy, i.e. temperatures that depend on the direction, the scale gmax(T) 
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dimensions 3 


££e — 2 trfjy + tr [D,-,Ao]" + m\ trAf) + — Ag(trAo ) 2 + ... 


(2.1) 


In lowest order one has a dimensionful coupling g 2 E = g 2 T + 0(g 4 ) and [38] 


m 2 E 


= (l + N f /6)g 2 T 2 +g 2 ^^ + 0(g 4 ), 


Xe = 9 ^ T + -- (2 ' 2) 


The dominant contributions to the thermodynamic pressure comes from the hard modes, and 
these contributions are completely perturbative, 


[jhard 


= r 4 (ci +c 2 g 2 +c 3 g 4 +c 4 g 6 + ...), 


(2.3) 


though the coefficients c >3 depend on the cutoff (A e) required to separate hard from soft scales. 
The soft contributions, on the other hand, are determined by the effective three-dimensional theory, 
and to three-loop order the result only involves the parameters gE and m E [ 0 ] 


Psoft/T = 


3n 


-m P - 


%TZ 2 


4 In 


A E 
2m E 


+ 3 I g\ m E ~ 


9 (89 11 


- 


1 


8 7T 3 V 24 6 


T ln2 + -7T" gE m E + 


6 


(2.4) 


The complete three-loop pressure of QCD is obtained by adding P = P hard + P s,>i \ and to achieve 
the maximal perturbative accuracy P hartl as well as the effective field theory parameters gE and m E 
are required to order g 4 . 

As long as T ^$> niE, this program can be extended to finite quark chemical potential. Besides 
modifying the parameters of the effective theory, there are also new, C-odd terms in the effective 
Lagrangian. The one with smallest dimension in nonabelian theories reads [39, ira 


‘ Sf * l)= *37c2^- /i * trA °. 


(2.5) 


In general the effects of these additional C-odd terms are small compared to the C-even operators. 


One quantity which is determined to leading order by the operator ([23]) is the flavour off- 
diagonal quark number susceptibility at zero chemical potential 


Xij - 


d 2 P 

dp-idp,;' 


( 2 . 6 ) 


When quark masses are negligible, all off-diagonal components are equal at /r,- = 0. Denoting 
them by the leading-order term involves a logarithmic term coming from the exchange of three 
electrostatic gluons and is given by 


(A 2 -i)(A 2 -4) m6 r2in i 

10 "' W g 


384A 


(2.7) 


is also the scale of magnetic instabilities @ 00 111- 

3 The relevant effective theory for time-dependent problems is still 3+1-dimensional and is given to leading order by 

the so-called hard-thermal-loop (HTL) effective action pi (35l |36|, |37|]. 
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Figure 2: The difference A P = P(T.p) P(T,()) divided by T 4 using the unexpanded three-loop result 
from dimensional reduction of Ref. |^8|] for fi/T = 0.2,..., 1.0 (bottom to top). Shaded areas correspond 
to a variation of it ms around the FAC-m choice by a factor of 2; dashed and dash-dotted lines correspond 
to the FAC-g and FAC-m results, respectively. Also included are the recent lattice data of Ref. | [f^ | (not yet 
continuum-extrapolated!) assuming Tf 0 = 0.49Aqcd- (From [|54|] ) 


where N is the number of colours. This vanishes in SU(2) gauge theory, but not in QCD, and also 
not in QED, where (in the ultrarelativistic limit) [ j42j ] 

e 6 2 1 

~- tT 2 In 

QED 


X 


( 2 . 8 ) 


24n 6 e' 

There has been for some time a discrepancy of the result ( |2.7[ ) with lattice results on off- 
diagonal quark-number susceptibilities. The authors of Refs. [ f43| . [44| | have obtained results in the 
deconfined phase that were far below those predicted by perturbation theory, and have interpreted 
this as new evidence for nonperturbative physics and the failure of weak-coupling methods. Most 
recent lattice results [45 have disproved the previous ones, and there is now agreement 

with the perturbative estimate at T > 2 T c . 

The complete three-loop pressure of QCD at finite quark chemical potential was recently cal¬ 
culated by Vuorinen [^8j]. Like the result at zero quark chemical potential, there is poor conver¬ 
gence and large renormalization scale dependence at realistic couplings, but again the apparent 
convergence can be improved importantly by applying the prescription of Ref. fll8 ] which keeps 
effective-field-theory parameters unexpanded. Figure || shows the result for AP = P(T.p ) — P(T.0) 
for Nf = 2 at several values of \i/T for which there are recent lattice data [4^]. The shaded re¬ 
gions correspond to variations of jU MS by a factor of 2 around a FAC value, and the dashed lines 
correspond to (two variants of) FAC jl ms- At T/To = 2, the highest value considered in [49j|. the 
FAC results exceed the not-yet-continuum-extrapolated lattice data consistently by ~ 10%, which 
is roughly the expected discretization error [[50]. When normalized to the free value APq instead 
of T 4 , the results would be essentially ;ti-independent and thus also very si mi lar to the Nf = 2 + 1 
lattice results of Ref. [|5~T|] as well as the quasiparticle model results of Refs. [|52j, 53]. 


3. Larg e-Nf limit of QCD and QED 

While the comparisons of weak-coupling results with lattice data are in remarkably good shape 
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Figure 3: Diagrams contribution to the thermodynamic potential of large-Ay QCD 


at temperatures a few times the deconfinement temperature, it seems desirable to have a cleaner 
testing ground for the resummation procedures required in applying weak-coupling techniques at 
realistic temperatures. Guy Moore [ ]55| ] has proposed to use the large-Ay limit of QCD and QED 
for this purpose. Large-Ay QCD is much simpler than large-A c QCD. In the limit Ay —> oo, A c ~ 1, 
g 2 Nf ~ 1, the relevant diagrams are those displayed in Fig. They involve a dressed gluon prop¬ 
agator which contains typical gauge-theory phenomena such as Debye screening for electrostatic 
modes, unscreened magnetostatic modes, complicated dispersion laws, Landau damping, and also 
plasmon damping. This is therefore a much richer theory than the large-A scalar field theories that 
are frequently used for testing resummations of thermal perturbation theory. 

In terms of the polarization tensor IL IV ' = I Iy^ + II^ t , where the matter part n^ at can be 
decomposed in a transverse and a spatially longitudinal piece, 11 / and IL , respectively, the thermal 
pressure reads 


„ TLjjtj 'In 2 T 4 f 2 T 2 n 

p = AAy ( —-+ 


+ N e 


180 

d 3 q 


6 

1 d(0 


llTt 2 


fj^l -J- 2 {[ n b + l] Imln(q 2 -w 2 + n T + n vac ) 
Imln (q 2 - or + n vac ) | 


, f f i i i t i o ® +n vac i T , q“ to + n v 

+ \ [ n b + j | Imln------ InrIn - 


q 2 - CO 2 


q 2 - CO 2 


+ 0(N j r 1 ). (3.1) 


The effective coupling is hidden in fl/ ^ and given by 


„2 


<?eff — 


tN f 


, QCD, 


k e 2 N f , QED. 


with exact one-loop beta function: 

i i MmTm) 

§1M £n(p') 6n 2 


(3.2) 


(3.3) 


There is no asymptotic freedom. In fact, there is a Landau singularity at exponentially large scales 
Al = /Ims^ 5 ^ 6 ^ 671 8ets ^ MS \ implying this theory exists only as a cutoff theory with A Cuto ff < A r . 
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Figure 4: Exact result for the interaction pressure at finite temperature and zero chemical potential at large 
Nf as a function of gg ff (/tMS = trT), compared to strict perturbation theory | |vs| . p6[ |57| ]. The tiny red band 
appearing for large values of the coupling for the exact result shows the cutoff dependence from varying 
the upper numerical integration cutoff between 50% and 70% of the Landau pole A l- The results of strict 
perturbation theory are given through order (dotted line), (dashed), (dash-dotted), and (dash- 
dot-dotted) where the renormalization scale Uvis is varied between \kT (line pattern slightly compressed), 
kT , and 2tzT (line pattern slightly stretched). The line labelled “FAC-m” indicates the scale chosen by the 
prescription of fastest apparent convergence for which the curves of and coincide. (From |^8]|) 


However, for the purpose of testing thermodynamic results this is no problem as long as the effec¬ 
tive cutoff provided by temperature and/or chemical potential T,p <C A| . It does however lead to 
certain technical intricacies, since Acutoff needs to be implemented such that Euclidean invariances 
are respected in order not to produce spurious singularities [ }55| . 56]. 

Comparing the available perturbative results with the “exact” large-A/ result obtained numer¬ 
ically [55, 56] shows that strict perturbation theory has the usual problems of poor convergence 
and large p -dependences (Fig. 4). The exact result has a curious nonmonotonic behaviour as a 
function of and only the small deviations from the free pressure at small coupling seem to be 
unambiguously predicted by perturbation theory. 

Figure [| shows the result of an optimized g 6 result with the prescription of keeping the 
effective-field-theory parameter me unexpanded. (The coefficient of the g 6 term has been extracted 
numerically with a 10% error shown by shaded bands.) The result agrees exceedingly well with the 
exact result even for the largest coupling where the ambiguity introduced by the Landau singularity 
is still under control. 

Perhaps even more remarkably, the nontrivial behaviour of the thermodynamic potential has 
recently been reproduced with comparable accuracy by using a 2-loop <t>-dcrivablc approximation 
to the entropy (where the perturbative accuracy is actually only order g 3 ) [59]. 

In Ref. [f57j] two of us have also evaluated the large-A/ pressure at finite chemical potential, 
which is displayed in the 3-d plot of Fig. |6| 

In Fig. |7j this result is compared with the perturbative result to order g 5 for two optimized renor¬ 
malization scales. While this shows a fairly large region where the perturbative results obtained 
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Figure 5: Exact result for the interaction pressure of large-A/ QCD at zero chemical potential [ |55| , |56| ] as a 
function of g^-fuvis = nT) or, alternatively, log !0 (Ai7'7r7'). The purple dashed line is the perturbative result 
when the latter is evaluated with renormalization scale /2ms = £fac = ne x ^~~^T ; the blue dash-dotted lines 
include the numerically determined coefficient to order (with its estimated error) at the same renormal¬ 
ization scale. The result marked “g^ rr = corresponds to choosing /2 ms such that the order-g® ff coefficient 
vanishes and retaining all higher-order terms contained in the plasmon term m 


90° 


N g (n 2 T 2 +in 2 ) 2 


- 0.001 



Figure 6: Exact result for the large-A/ interaction pressure P — Pq normalized to N s {k 2 T 2 + p 2 ) 2 as a 
function of gg ff (/2Ms) with /2 ms = k 2 T 2 + p 2 , which is the radial coordinate, and (p = arctan(7iT / p). 
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Figure 7: Percentage errors of the perturbative result for the interaction part of the pressure to order y^ |f 
in the larg e-Nf limit for two choices of Pms : Fastest apparent convergence of P as well as m\ (FAC-m), 
and of g/ : (FAC-g). The brightest area corresponds to an error of less than 1%, the darkest ones to an error 
of over 100%. The ratio of chemical potential to temperature increases from top to bottom according to 
(p = arctan(7T T/Li) so that 90° corresponds to high temperature and zero chemical potential, and 0° to zero 
temperature and high chemical potential. The coupling is given in terms of gg ff (/tMs) at Ams = sjtt 2 T 2 + p 2 . 
(From [Q) 


through dimensional reduction is a good approximation, it also shows that there is a breakdown of 
dimensional reduction at very small temperature T < O.ljU. In this regime it becomes necessary 
to evaluate all infrared sensitive diagrams such as ring diagrams in four dimensions. Combining 
such procedures with the existing analytical results to three-loop order should be able to fill in the 
existing gap of perturbative results on the pressure for very small temperatures and high chemical 
potential [|60l|. 

A three-loop result for the QCD pressure does exist at exactly zero temperature. It is due to 
Freedman and McLerran [67j. 52, 63] and it has recently been updated by Vuorinen [48]. In the 
larg e-Nf limit it reads 


P-Po 

Seff 

1 § e ff ^ - AMS n ocoo 

stn 

N g p* 

t= o 32 tt 4 

411 r\ 111 v/» J • • • 

2n 2 3 ju 

1287T 6 


+ ^(grfflngeff)- (3-4) 
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Figure 8: The interaction part of the pressure at zero temperature and finite chemical potential as a function 
of tQt^AMS = M) or, alternatively, log| (l (A| //i), compared with the perturbative result of Freedman and 
McLerran | fi7j |63| [ to order and our numerically extracted order-g® ff result, both with renormalization 
scale in the perturbative results varied around ft ms = /r by a factor of 2. The coloured bands of the g^-results 
cover the error of the numerically extracted perturbative coefficients. (From [|57|l) 


The exact (but numerical) large-A/ result of Ref. [ p7| | in fact provides a check of these results 
as well as a numerical determination of the next terms to order and ln(y e ff): 


P-Po 


geff> 7 ’=0,MMS=/r 


2n 2 

[3.18(5)ln-5- + 3.4(3)] 

geff 


geff 

2048 tt 8 ' 


(3.5) 


Figure |8] compares the exact and the perturbative results for the large -Nf pressure at zero 
temperature and high chemical potential, with the perturbative results evaluated for three values of 
Mms/M- 


4. Non-Fermi-liquid behavior 


From the leading-order interaction term in the thermal pressure [64] 


P-Po = ~ N g 



g|f , 

32 


(4.1) 


one would expect that the entropy density 5? = at small T <C. f should start as 

y-S'o = -N g ^f 2 T + ... (4.2) 

Flowever, for sufficiently small T, namely when T <C gjJL, the interaction part of the entropy 
density extracted from the exact large-A/ result (Fig. |9|) has even a different sign. The contributions 
to ( |3.ip which do not involve the Bose distribution function n/, do indeed behave in accordance to 
the result ( |4.1[ ), but those with the function /i/ ; that one would expect to be negligible at very small 
temperature turn out to become the dominant ones, and they contribute terms which are nonlinear 
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Figure 9: The interaction part of the entropy at small T/jl for = ju) = 1, 4, and 9. The “non-/;/,” 

contributions (dash-dotted lines) are negative and approximately linear in T with a coefficient agreeing with 
the exchange term in the pressure at small coupling; the “ rq, ” contributions (dashed lines), which are 
dominated by transverse gauge boson modes, are positive and nonlinear in T such that the total entropy 
exceeds the free-theory value at sufficiently small T /p. 


in T, implying a strong deviation from the usual behaviour of Fermi liquids [|65p. This non-Fermi- 
liquid behaviour is caused by long-range interactions, namely weakly screened magnetic modes. 
At very small frequencies, the latter have a dynamical screening length K = [nm^to/A] 1 ■ and 
it has been found long ago by Holstein, Norton, and Pincus [ j66{ ] that this manifests itself in the 
appearance of an anomalous contribution to the low-temperature limit of entropy and specific heat 
proportional to cxT\nT '. 

While this effect is perhaps too small for experimental detection in nonrelativistic situations, 
it drew renewed theoretical attention more recently [ 


68 , 


after the detection of non-Fermi- 
liquid behavior in the normal state of high-temperature superconductors [f70j] and in other systems 
of strongly correlated electrons, which may be due to effective gauge field dynamics (see also 

00 ®)- 

In deconfined degenerate quark matter, the analogous effect can more easily be important 
because the larger coupling constant a s together with the relatively large number of gauge bosons 
increases the numerical value of the effect by orders of magnitude. In contrast to the case of a high- 
temperature quark-gluon plasma, chromomagnetostatic fields are expected to remain unscreened 
in the low-temperature limit [21] and thus lead to the same singularities in the fermion self-energy 
that are responsible for the breakdown of the Fermi-liquid description in the nonrelativistic electron 
gas considered in [j66j]. 

An important consequence of such non-Fermi-liquid behavior in quantum chromodynamics 
(QCD) is a reduction of the magnitude of the gap in color superconductors i fTij , [25], [26], 27 ] which 
on the basis of weak-coupling calculations are estimated to have a critical temperature in the range 
between 6 and 60 MeV [74j[. Quark matter above this temperature has long-range chromomagnetic 
interactions that should lead to an anomalous specific heat with possible relevance for the cooling 

]. However, in Ref. [ |76| | 


of young neutron stars as pointed out by Boyanovsky and de Vega [75 
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these authors claimed that the aT\nT 1 term in the specific heat as reported in [66, 68, |69|] would 
not exist, neither in QCD nor in QED. Instead they obtained a aT 3 In T correction to the leading 
ideal-gas behavior, which by renormalization-group arguments was resummed into a 7 3 10<a> cor¬ 
rection as the leading non-Fermi-liquid effect on the specific heat. 4 At low temperatures, such a 
contribution would be rather negligible compared to standard perturbative corrections to the ideal- 
gas result oc T. 


In Refs. \fn\ [78Q we have recently been able to resolve these contradictory results and shown 
that the aT\nT 1 term as obtained in []66j, [SSl |69|] is indeed correct. In addition to the coefficient 


of this leading logarithm, we have calculated the coefficient under the log as well as higher terms 
in the low-7 series, which come with fractional powers of 7, 


N, 


sln^ 2 T 

36n 2 

2 2/3 r(|)C(f 


1*-. r“) | n 6 (r^\ 

ln ^ F - 2 + 7£- F ?( 2 ) 


3 '"’' (ge«r) i,3 T^ + 8 i) 2 ‘Cr('°)C ( '°) fa ^) 2 / 3 7 . 7/ 3 


+- 


9V/37T 11 / 3 
2048-2567T 2 -367r 4 + 37r 6 
5407T 2 


27\/37r 13 / 3 
l n M^_4.3493485... 


+ 0(r n / 3 ). (4.3) 


This expansion is systematic for 7 <C y e irF- but breaks down at 7 ~ ^ e n/J; for 7 A> g e ff/i it 
in fact has to switch to the perturbative series that can be obtained from dimensional reduction. 
A result which covers all 7 <C jU can be obtained by resumming the complete nonlocal HDL self 


energies and is given by the compact expression [78] 


—{y-y°) = - 


N t 


247T 2 


1 /■” , dn b {q 0 ) r , • 

lie Jo dq "-3T-L dqq 


2Ira in (TzilEEl) +Im 

+0(gi^ 2 T). 


12 \ 2 2 

T~% J \ T~% 


(4.4) 


For 7 <C gft this resums the above low-7 series which is dominated by magnetic effects and for 
g c(( jl <C 7 <C 7t connects smoothly to the dimensional reduction result 


iy g 




Sln^ 2 T g e V 3 

8 ?r 2 12 tt 4 


(4.5) 


where the 0(g 3 ) contribution is from electrical Debye screening. The numerical evaluation of ( fOj ) 
together with the first few orders of the low-temperature series and the perturbative result from 
dimensional reduction is shown in Fig. 11- 

Figure [TT] displays the resulting specific heat normalized to its ideal-gas value. The lines 
marked “QED” correspond to g e ff = 0.303 or (Xqed ~ 1/137, and the results for y e rr = 2,3 corre¬ 
spond to a s ~ 0.32,0.72 in two-flavor QCD. (Recall that g 2 ff = g 2 Nf/ 2.) While in QED the effect 
is tiny (the deviations from the ideal-gas value have been enlarged by a factor of 20 in Fig. [TT] to 
make them more visible), in QCD we find that there is an interesting range of T/\i where there 

4 Resummation of the a71n7 _1 term along the lines of Ref. |7^] would have led to a 7’ 1 +°( o 0 term instead. 
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Figure 10: The function 6(T /(g e ffM)) which determines the leading-order interaction contribution to the 
low-temperature entropy. The normalization is such that © = — 1 corresponds to the result of ordinary 
perturbation theory. The various dash-dotted lines give the first few orders of the low-temperature series for 
the entropy; the perturbative result from dimensional reduction to order g 2 ff and is given by the full and 
the long-dashed curves, resp. (from [[78]]) 


is a significant excess of the specific heat over its ideal-gas value, whereas ordinary perturbation 
theory [0 ] would have resulted in a low-temperature limit of ( S\f ( S' {3 = 1 —2a s /n. At sufficiently 
low temperature, may even become much larger than 1. As discussed in Ref. [|79|], the large 

logarithm that appears in this limit is actually stable against higher-order corrections. However, 
in QCD (though not in QED) the growth of is limited by the appearance of a supercon¬ 

ducting phase at sufficiently small T. According to Ref. [74], the critical temperature for the color 
superconducting phase transition may be anywhere between 6 and 60 MeV, so with e.g. a quark 
chemical potential of jU = 500 MeV the range T/fJ, > 0.012 in Fig. [IT] might correspond to normal 
quark matter with anomalous specific heat. While the effect remains small in QED, it seems con¬ 
ceivable that the anomalous terms in the specific heat play a noticeable role in the thermodynamics 
of proto-neutron stars, in particular its cooling behavior in its earliest stages before entering color 
superconductivity [8(ij, |8T|]. 

The cooling of (proto-)neutron stars with a normal quark matter component is not only sen¬ 
sitive to non-Fermi-liquid effects in the specific heat. The latter are in fact overcompensated by 
two powers of a v ln (m/T) appearing in the neutrino emissitivity [ j82| ] which come from logarith¬ 


mic singularities in the group velocities of quark quasiparticles. While Ref. [82] assumed a single 
scale for the two kinds of non-Fermi-liquid logarithms, a recent explicit calculation [ 83 ] showed 
that these can be very different. 

The scale of the logarithm in the specific heat is approximately given by log(0.282m/r) with 
m 2 = Nf(gp /In) 2 . At zero temperature, the group velocity of quark quasiparticles near the (would- 
be) Fermi surface (£ = E — ft) reads [83] 




.07 m 


0((e/m) 2 / 3 ) 


(4.6) 
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Figure 11: The HDL-resummed result for the specific heat normalized to the ideal-gas value for g e ff = 2 
and 3 corresponding to a s ~ 0.32 and 0.72 in two-flavor QCD, and g e ff ~ 0.303 for QED. The deviation of 
the QED result from the ideal-gas value is enlarged by a factor of 20 to make it more visible. 


(v 1 - Wc» 



logto (e/m) 

Figure 12: The group velocity of quark quasiparticles plotted as (vT 1 — 1 )/{g 2 Cf) over log 10 (£’ — p/m) at 
zero temperature. 


which at small nonzero temperature T <C gjx is bounded by 

v g - 1(o) = l + 01"^ + o ( (r/ m ) 3 ). (4.7) 

Evidently the scale of the latter logarithm differs from that in the anomalous specific heat by more 
than a factor of 30. 

Higher terms in the small £ expansion of the group velocity v“ (e) also involve fractional 
powers, and a complete result for general £ can again be obtained by HTL resummation. The latter 
is shown in Fig. [H| together with the first few terms of the small-£ expansion. 
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5. Conclusions 

There has recently been substantial progress in applying weak-coupling methods to the ther¬ 
modynamics of deconfined QCD, and we have made a case that then - generally poor apparent 
convergence can be greatly improved by suitable resummations. Results obtained from the effec¬ 
tive field theory of dimensional reduction remain predictive for temperatures down to a few times 
the deconfinement temperature (rather than many orders of magnitude higher as previously con¬ 
cluded). This also holds true for nonzero chemical potential fi ~ T, and in fact comparisons with 
the available lattice data show reasonable agreement already at T ~ 27). 

For temperatures T ~ gfi or smaller, the perturbative expansion obtained in the dimensional 
reduction framework eventually breaks down because of non-Fermi-liquid effects. Those can be 
calculated analytically by HDL resummation and give rise to anomalous terms in the specific heat 
and related quantities. While tiny in QED, such effects may become important in the thermody¬ 
namics of neutron stars with a normal quark matter component. 
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